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I.   INTRODUCTION 

Fuel  tanks  comprise  the  greatest  portion  of  vulnerable  target 
volume  in  modern  military  aircraft.   This  degree  of  vulnerability  of  the 
fuel  tanks  is  of  paramount  importance  in  determining  aircraft  surviva- 
bility.  Fuel  tanks  which  are  subjected  to  ballistic  impact  can  sustain 
severe  damage  as  a  result  of  the  phenomenon  known  as  hydraulic  ram.   For 
our  purposes,  hydraulic  ram  is  defined  as  the  production  of  pressure 
loadings  by  ballistic  impact  and  their  effects  on  the  fuel  tank  and 
its  components.   The  term  hydraulic  ram  is  a  misnomer  in  that  it  sug- 
gests a  quasi- static  compression  process.  Actually  the  mechanism  of 
hydraulic  ram  consists  of  several  phases.   The  shock  phase  is  generated 
when  the  projectile  penetrates  the  tank  wall.   Projectile  penetration 
produces  a  stress  riser  in  the  tank  wall  from  which  cracks  can  propa- 
gate radially.   Shock  wave  formation  in  the  fuel  due  to  projectile 
penetration  into  the  fluid  causes  very  high  local  pressures  which  may 
be  sufficient  to  cause  failure  of  the  entry  wall  in  the  neighborhood  of 
the  entry  point.   As  the  projectile  moves  through  the  fluid  it  leaves  a 
vapor- filled  void  caused  by  flow  separation  from  the  projectile  surface. 
During  this  so-called  cavity  phase  the  projectile's  kinetic  energy  is 
transferred  to  the  fluid  by  cavity  formation.   Cavity  growth  is  subse- 
quently defeated  by  hydrostatic  pressure  and  the  cavity  collapses. 
The  collapse  is  not  total  because  of  fluid  vapor  compression.   The 
cavity  oscillates  until  equilibrium  is  reached  in  the  fluid.   Cavity 
oscillations  produce  corresponding  pressure  loadings  on  the  tank  walls 
which  can  cause  catastropic  failure  of  tank  components.   The  exit  wall 
of  the  tank  is  also  heavily  loaded  as  the  projectile  passes  through. 

The  purpose  of  the  investigation  to  date  has  been  to  develop  an 


understanding  of  the  hydraulic  ram  phenomenon.   Fuel  tank  response  to 
ballistic  penetration  cannot  be  fully  understood  without  comprehension 
of  the  individual  hydraulic  ram  components  and  their  interaction.   This 
understanding  can  be  obtained  by  conducting  experiments  and  developing 
analytical- numerical  tools  which  predict  the  events  that  occur  during 
the  damaging  processes.   With  this  understanding  and  the  analytical- 
numerical  tools  available,  fuel  tanks  which  minimize  hydraulic  ram  can 
be  designed.   Once  fuel  tanks  are  protected  against  this  type  of  threat, 
aircraft  vulnerability  will  be  lessened. 


II.   ANALYTICAL  AND  NUMERICAL  STUDIES 

A.   A  Study  of  the  Wall  Motion,  Treated  as  a  Single  Degree  of  Freedom 

System  (Development  of  Piston  Theory) 

Consider  a  typical  portion  of  a  fuel  tank  wall  treated  as  a  single 
degree  of  freedom  system  as  shown  in  Fig.  II-A-1.   The  area  of  the  por- 
tion of  the  wall  is  A  and  the  total  mass  is  M.   The  resistance  to 
displacement  provided  by  the  surrounding  wall  is  modeled  by  the  spring, 
whose  constant  is  K.   The  equation  governing  the  elastic  motion  of  the 
wall  is  assumed  to  be 

Mw  +  Kw  =  Ap(t)  ( II-A-1) 

where  p(t)  is  the  total  time  dependent  pressure  on  the  wall  due  to  the 
incident  pressure  waves  in  the  fluid  and  to  the  waves  reflected  and 
radiated  from  the  wall  and  w  is  the  wall  displacement.   The  dots  indi- 
cate derivatives  with  respect  to  time  t.   Assuming  that  the  fluid 
satisfies  the  linear  wave  equation  and  applying  a  step  pressure  pulse 
normal  to  the  wall  as  the  incident  wave  gives 

p(t)  =0    t  <  0 


p(t)  -  p0  -  ^   t  ,  o  ("-A-2> 


where  p  is  the  pressure  in  the  incident  wave,  p  is  the  fluid  density, 
and  f>   is  the  velocity  potential  of  the  reflected  and  radiated  waves. 
The  subscript  i  denotes  the  value  of  f)   at  the  wall.   The  potential  fi 
is  governed  by  the  wave  equation 

c2fi     =   J5  (II-A-3) 


where  c  is  the  acoustic  velocity  in  the  fluid  and  primes  denote  deriva- 
tives with  respect  to  x,  the  axial  space  coordinate  normal  to  the  wall, 
At  the  wall  the  velocity  of  the  fluid  is  equal  to  the  velocity  of  the 
wall .   Thus , 


w  =  v.  +6.  (II-A-1+) 

1   ri 


where  v  is  the  velocity  of  the  incident  wave  and  is  equal  to  p  /(pc). 
i  o 

Equations  (l),  (2),  (3),  and  (k)  govern  the  behavior  of  the  fluid- 
structure  interaction  process.  The  solution  to  Equation  (3)  for  f>  can 
be  given  in  the  form 

<h   _   i^(*+ct)/a   B  U(x-ct)/a 

P  ~  ae  +  Pe  (II-A-5) 

where  a   and  (3  are  arbitrary  constants,  X  is  a  scalar,  and  a  is  the 
value  of  x  at  the  wall.   The  first  quantity  on  the  r.h.s.  of  Equation  (5) 
represents  waves  traveling  away  from  the  wall,  and  the  second  quantity 
represents  waves  traveling  toward  the  wall.   Note  that  for  the  waves 
traveling  away  from  the  wall,  i.e.,  the  reflected  and  radiated  waves, 

t      t 

P  =  cfi  (II-A-6) 

Thus, 

w  =  v.  +  -  0 
l   c  r 

according  to  Equation  (k) .      Hence 

p  =  c(w-v.  )  (II-A-7) 

Substituting  fi   given  by  Equation  (7)  into  Equation  (2)  gives 


P(t)  =  PQ  +  pc(vi-w)     t  >  0 


(II-A-8) 


Examination  of  Equations  (l)  and  (8)  reveals  that  the  variables  P   and 
w  have  now  been  separated,  i.e.,  the  response  of  the  wall  can  be  compu- 
ted based  upon  the  pressure  and  velocity  of  the  incident  pressure  wave. 
Application  of  Equation  (8)  to  the  more  general  two-dimensional  wall 
response  is  known  as  the  piston  theory  (Refs.  1  and  2).   Substituting 
Equation  (8)  into  Equation  (l)  leads  to 


•or'  2    P° 

w  +  2Cw  +  uo  w  =  — 

*  m 


(II-A-9) 


where 


2   K 

U)  -   — 

M 


b   2muu 


p  =  p  +  pcv.  =  2p 
*o   -^o      1     o 


M 
m  =  A 


The  solution  to  Equation  (9)  can  be  given  in  the  form 


w  =  w  ,  J  1-e 
st  ^ 


^c+v/r-i  p^y^it  _  q-Jq2-i  0-*Jf- 


■*-^.   e 


U/ 


c£-i 


$ 


/C2-l 


(II- A- 10a) 


w  =  w  ,  e  *   I  e  ' 

st 


2^T 


J^l     t  _  e-wj£-l 


(II- A- 10b) 


.ypTZ  t  _  e-w/?I7 


where 


w 


Ap"  p 

yo         o 

st   "     K  2 

muo 


(11-A-lOc) 


when  the  wall  is  at  rest  prior  to  the  application  of  the  step  pressure. 

If  the  contribution  of  the  spring  is  neglected,  we  have  a  one- 
dimensional  pressure  wave  reflecting  off  of  an  infinite  wall.   Hence, 

1  PC 
oi  =  0,  £  =  00>  C00  =  o  —  •   This  problem  was  considered  in  Ref.  3,  where 

the  equation  for  the  pressure  on  the  wall  was  taken  as 

pet 

(Il-A-I0d) 


p(t)  =  pQ  e 


m 


Typical  values  for  c,  p  and  m  for  hydraulic  ram  in  fuel  tanks  are 


c  =  53>600  in. /sec 


p  =  7.64  x  10"5  #-sec  /in 


2/.   If 


-5     2    3 
m  =  2.59  x  10   #sec  /in.    (0.1  in.  aluminum  wall) 


=  3.64  x  10"  #sec2/in.3  (0.5  in.  steel  wall) 


Therefore 


■£2  =  0.158  x  10  —  (0.1  in.  aluminum  wall) 
m  sec 


-^  =  0.01125  x  10  —  (0.5  in.  steel  wall) 
m  sec 


Consequently,  when  u>  <  15000  rad/sec  for  the  0.1  in.  aluminum  wall 
and  uo  <  1000  rad/sec  for  the  0.5  in.  steel  wall,  £  >  >  1,  i.e.,  the 


system  is  heavily  damped,    and  Equations   (10a  -    10c)    simplify  to 


uut  U)t 

,  uu   \   [   "   2Q         2£     -2£uut| 
Wst     l2CUe 


(II- A- lib) 


orb         uut 
2£         2£     -2£uut' 
Po  \X  ~   e  +  e        e  j 

A  plot  of  —  versus   time   is  given  in  Fig.    II-A-2a  for  —  =  0.2,   0.1, 
p  m 

6     1 

0.05,    and  0.01  x  10     ,    and  1  |isec.      In  this   figure   the   assumption  is 

o        sec 

/:   -. 

made  that  —  <  100,  i.e.,  uu  <  1,000  rad/sec  for  -££  =  0.01  x  10  — 
pc  '        m  sec 

and  oo  <  U.U70  for  —  =  0.2  x  10  .   Under  these  conditions,  Equation 

m  sec 

(lie)  for  p  simplifies  to  Equation  (lOd),  provided  t  <  500  u-sec.   Curves 
of  this  type  are  also  given  in  Ref.  3,  Fig.  8,  for  a  1  inch  plexiglas 
wall  and  a  0.063  inch  aluminum  wall  with  an  exponentially  decaying  inci- 
dent pressure  wave. 

For  t  >  500  (isec  the  elastic  property  of  the  wall  begins  to  become 

effective  and  reduces  the  wall  velocity.   This  causes  the  pressure  to 

2 
build  up.  A  plot  of  |-  is  given  in  Fig.  II- A- 2b  for  —  =  100,  20,  10 

1  °  P° 

and  5  >  and  100  u^sec  £  t  £  100,000  M,sec.   For  these  conditions, 

SGC 

Equations   (lla  -   lie)    simplify  to 

uut  /uut 


Wst    l1"6     "*]     >W  =  Wstl2Tje     ^      ,   P  =P0  U-e2'  e"' 


(11-A-12) 

Examination  of  Figs.  II-A-2a  and  2b  reveals  that  for  the  values  of 
the  parameters  considered,  the  inertial  effects  and  the  elastic  effects 
occur  in  essentially  separate  periods  of  time.   This  is  due  to  the 

8 


assumption  on  uj.  When  the  stress  wave  strikes  the  wall,  the  wall  de- 
flects as  if  there  were  no  elastic  effects.   After  approximately  500- 
1000  M-sec  the  wall  starts  to  slow  down  due  to  the  elastic  effects,  the 
pressure  builds  up,  and  eventually  the  wall  comes  to  rest.   The  motion 
of  the  wall  is  that  of  a  heavily  damped  oscillator*.   This  phenomenon 
holds  provided 

™L  <  100  -A- 

m       sec 


For 


■££  =  0.01  x  10     —      ,      oj  <  1000  rad/sec 
m  sec      7  ' 


The   lowest  natural   frequency  of  a  square  plate   of  width  L  is 

2 


2tt 

0)    = 


V12(l-^ 


L2     Vl2(l-u2)p  (II- A- 13) 


where  E  is  Young's  modulus,  h  is  the  wall  thickness,  and  u  is  Poisson's 
ratio.   For  both  aluminum  and  steel  walls 


uu  ^  1.2  x  105  f^\   rad/sec 


L2i  (II- A- Ik) 


For 

L  =  20  in.   and  h  =  0.5  in.     ,     u)  «  15 0  rad/sec 

For 

L  =  10  in.      and     h  =  0.5   in.  ,  uu  sa  600  rad/sec 


*The  natural  frequency  for  critical  damping  is  79,000  rad/sec   (12,600 
Hz)    for  the  0.1  aluminum  wall,    and  is   5,620  rad/sec  (884  Hz)    for  the 
0.5   in.    steel  wall. 


Thus,  a)  <  1000  rad/sec  is  a  realistic  limit  for  the  fundamental  frequency 
of  fuel  tank  walls. 

Since  the  elastic  response  takes  place  relatively  late  in  the  inter- 
action process,  the  question  arises  as  to  the  possible  validity  of 
treating  the  stress  wave  as  an  impulse  loading.   For  example,  suppose 
the  incident  pressure  pulse  at  the  wall  is  of  the  type  shown  below 


P  >N 


■>   t 


The  solution  for  the  wall  displacement  (provided  C  >   l)  for  t  >  t  is 


w  = 


■il  +' 


K 


2  Cm 


(t-T)uu 
2C 


t    >    T 


w  w 

T  T 

+ 


LkC 


2        2Cw 


Uu(t-T) 

2C 

e  * 


-2£w(t-T) 


(II-A-15) 


where  w  and  w  are  the  displacement  and  velocity  of  the  wall  at  t  =  t. 
The  first  term  in  Equation  (15)  is  the  long  time  response;  the  second 
term  damps  out  very  quickly.   If  the  stress  wave  pulse  is  treated  as  an 
impulsive  loading,  the  displacement  and  velocity  at  t  =  t  are 


w  =  0 

T 


W   =  T 

t   m 


Hence 


w 


V 


impulse    2  Cum 


u>(  t-  t  ) 

2C 


-  e 


Oj(t-T) 

2C 


-2Cu)(t-T) 
e         y  (II- A- 16) 
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The   correct   solution  for  w  using  Equations   (10a)    and  (10b)    to  define 


w     and  w     is 

T  T 


UOT 

w     =  — «     1-e 


t  2  \  /      2£uu 

moo 


UOT  UOT 

Po  "   21  2C 


and 


„        .   -\     -  %ll  ]_      /      |I     -2£U»T\     i%H      -2C.(t-T 


w  f«  — 5      <  1-e  e  to ^   ll-e 

moo2       1\  J  h£ 


t    >     T 


provided   Q  >  >  1.      When  t  is   such  that 

UUT 

"    2£  OJT 

e  ~  i  -   2Q 


i.e.,    t  <  <  — p    ,   Equation  (17)    simplifies  to 

moo 


(II- A- 17) 


(yj(t-T)  _,.    U0(t-T)    _r  /.   A 
^2C       1   "  ^T    2f    -2£w(t-T) 
e     "    +  — ^  e      e   *   e        / 


Po        2£  2£  t  > 


2C^m  y  kg 


(II-A-18) 


Thus,  comparing  Equations  (l6)  and  (l8)  reveals  that  the  long  time 
response  is  correctly  approximated  by  the  impulse  assumption,  but  the 
short  time  response  is  inaccurate. 

The  conclusions  that  can  be  drawn  from  this  analysis  are: 
1.   The  wall  responds  to  the  pressure  wave  as  a  very  heavily 

damped  oscillator  provided  the  natural  frequency  of  the  wall 
is  less  than  1000  rad/sec  (rJ.50  Hz). 


11 


2.  The  pressure  on  the  wall  drops  off  to  a  small  percentage  of 
the  incident  pressure  within  the  first  100  [isec.      It  begins  to 
build  up  when  t  >  1000  u.sec.   Thus,  the  reflected  pressure  from 
the  wall  is  essentially  the  same  as  the  reflected  pressure  from 
a  free  surface  for  the  period  100  \isec   <  t  <  1000  t-isec. 

3.  If  the  incident  pressure  is  on  for  a  finite  length  of  time  t 

such  that  t  <  <  —  (  —  =  0.158  x  10  for  a  0.1  in.  aluminum 

pc  \m  sec 

wall,  -££  =  0.01125  x  10  —  for  a  0.5  in.  steel  wall)  the  long 
m  sec  / 

time  response  is  accurately  predicted  but  the  short  time  response 
is  inaccurate. 
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Figure  II-A-1   Idealization  for  Single  Degree  of  Freedom  Analysis 
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B.  Entry  Wall  Response  to  the  Shock  Phase 

A  preliminary  study  of  the  structural  response  of  the  entrance  wall 
of  a  fuel  tank  during  the  initial  shock  phase  of  the  hydraulic  ram  has 
been  conducted.   The  solution  for  the  pressure  in  the  fluid  derived  by 
Yurkovitch  (Ref.  k)   was  reduced  to  a  simple  analytical  form  for  use  in 
the  structural  analysis  digital  computer  program  SATANS*.   The  form 
selected  for  the  pressure  p  at  a  radius  r  and  time  t  was 


P(r,t)  =  p^-jVT-;  [Hhr]      o  <  t  <  t1c1 


(H-B-la) 


(H-B-lb) 

where  1  ^  {  }  >  0.   This  form  corresponds  to  a  peak  pressure  that  varies 

t  "m 
as  p*  —    with  a  linear  decay  behind  the  front  to  zero  pressure.   The 


W- 


location  of  the  front  is  given  by  Cp  [  -^—     .   For  a  0.50  caliber  project- 
ile impacting  Hp0  (Fig.  22  of  Yurkovich's  report) 

p*  =  150  ksi 

t  =  tp  =  t„  =  10  jasec 

C  =  1.22 

C2  =  1.7  in. 

m  =  1.2 


*SATANS  computes  the  transient,  geometrically  nonlinear  response  of 
elastic,  axisymmetric  plates  and  shells  subjected  to  arbitrary  loads 
It  is  described  in  detail  in  Refs.  5  and  6. 
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n  =  O.87 

q  =  0.68 
The  pressure  loading  given  by  this  model  is  shown  in  Figure  II-B-1  and 
was  applied  to  a  circular  aluminum  plate  with  a  radius  of  10  in.  and 
a  thickness  of  0.1  in.   The  outer  edge  of  the  plate  was  assumed  to  be 
clamped.   The  results  for  the  plate  displacements  computed  by  SATANS 
are  shown  in  Fig.  II-B-2  for  t  =  5,  10,  20,  30,  and  50  [asec.   The 
maximum  membrane  force  in  the  plate  was  approximately  60,000  lb/ in. 
and  occurred  at  the  center  of  the  plate  at  t  «  35  u.sec.   The  maximum 
stress  caused  by  the  force  and  by  the  bending  moment  was  600,000  psi 
and  720,000  psi  respectively. 

The  analysis  for  the  response  of  the  entrance  wall  during  the  shock 
phase  discussed  above  neglects  several  important  features.   First,  it 
does  not  take  into  consideration  the  fact  that  the  motion  of  the  wall 
relieves,  to  some  extent,  the  pressure  acting  on  the  wall.   Second, 
the  very  large  stresses  developed  in  the  wall  are  considerably  higher 
than  the  yield  stress  and  the  rupture  stress  of  aluminum.  Hence, 
yielding  and  cracking  will  occur.   Third,  the  effects  of  shear  defor- 
mations and  rotary  inertia  have  not  been  considered.   Of  these  three 
features,  the  effect  of  wall  motion  upon  the  loading  pressure  is  consi- 
dered to  be  of  prime  importance  in  the  initial  phase  of  this  program. 
An  exact  method  for  determining  the  total  fluid  pressure  within  a 
closed,  flexible  tank  has  been  developed  by  Lundstrom  (Ref.  6).   This 
model  was  considered  as  the  algorithm  for  the  pressure  loading. 
However,  due  to  difficulties  encountered  in  the  programming  of  the 
model,  Lunstrom  recommended  that  the  piston  theory  be  used  in  place 
of  his  exact  model.   Piston  theory  is  described  in  II-A.   Application 
of  this  theory  uncouples  the  structure  and  fluid  equations  and  allows 
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the  damped  response  of  the  wall  to  be  computed  using  the  incident 
pressure  as  the  load  on  the  wall.   This  approximation  has  been  in  use 
since  the  early  19^0' s  when  it  was  applied  to  the  study  of  the  effect 
of  underwater  explosions  on  ship's  plates. 

The  theory  is  based  upon  the  assumption  that  the  normal  pressure  on 
the  wall  is  given  by  Equation  II-A-8,  i.e., 

p(t)  =  p  +  p  (v.  -  w) 
o    c   1 

which  is  the  equation  for  one- dimensional  fluid  flow.   Since  SATANS  did 
not  contain  the  damping  term  pew,  it  was  necessary  to  modify  the  code 
to  incorporate  this  procedure. 

The  results  for  the  response  of  the  same  circular  plate  subjected 
to  the  same  pressure  pulse  described  above  but  now  with  piston  theory 
added,  are  presented  in  Fig.  II- B- 3  at  several  times.   Comparing  the 
two  figures,  II-B-2  and  3,  reveals  that  accounting  for  fluid- wall  inter- 
action has  a  considerable  effect  on  the  wall  response.   Nevertheless, 
for  this  example  the  wall  stresses  are  still  above  the  yield,  i.e., 
the  maximum  membrane  force  is  8,900  lb/in.  and  occurs  at  the  center  at 
t  =  15  M-sec,  the  maximum  bending  moment  is  V78  in. -lb/in.  and  occurs  at 
r  =  1  in.  at  t  =  15  psec,  the  maximum  stress  is  327  ksi  and  occurs  at 
r  =  1  in.  at  t  =  15  u-sec. 
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FIGURE   H-B-2      NORMAL    DISPLACEMENT      vs.      RADIAL 

DISTANCE    WITHOUT    FLUID    INTERACTION 


20 


E 
a> 
o 
o 

Q. 
(A 


I  2  3 

Radial     Distance  (in.) 


FIGURE     I-B-3       DISPLACEMENT    OF    FUEL    TANK 

WALL,  WITH    PISTON     THEORY 


21 


C.   Modification  of  the  Northrop  Internal  Blast  Code  (BR-l)  to  Account 
for  Hydraulic  Ram  Loading  and  Fluid- Structure  Interaction 
The  Northrop  Corporation  is  developing  a  finite  element  digital 
computer  code  for  predicting  the  transient,  inelastic,  large  deflection 
response  of  flight  vehicle  structures  subjected  to  internal  blast  loading 
and  fragment  penetration.   Ball  has  been  following  the  development  of 
this  code  with  the  objective  of  modifying  it  to  account  for  the  hydraulic 
ram  loading  and  fluid- structure  interaction  that  occurs  when  bullets 
and  metal  fragments  penetrate  aircraft  fuel  tank  walls.   The  main  features 
of  the  program  are  summarized  below: 

1.  The  elements  are  the  isotropic  flat  rectangular  plate  and  the 
straight  beam.   The  cross- section  of  the  beam  element  is  either  rectan- 
gular, I,  or  T-shaped. 

2.  The  Bernoulli- Euler  hypothesis  is  used  for  both  the  plate  and 
the  beam  element. 

3.  The  nodes  of  the  plate  element  are  at  the  four  corners.   The 

degrees  of  freedom  at  each  node  are  the  three  displacements  u,  v,  and  w 

and  two  rotations  9  and  9  .   A  linear  edge  displacement  state  is  assumed 

x      y 

for  the  inplane  displacements  u  and  v,  and  a  cubic  function  is  assumed 
for  the  normal  displacement  w. 

k.      Large  deflection  strain-displacement  relationships  are  used  for 
the  plate  element  that  are  valid  for  strains  as  large  as  10$  and  rota- 
tions as  large  as  20  . 

5.  Gaussian  quadrature  is  used  for  the  integration  process.  Three 
points  in  each  direction  are  used. 

6.  The  displacement  state  assumed  for  the  beam  element  leads  to  an 
axial  force,  two  transverse  shearing  forces,  and  a  uniform  torsion  over 
the  length  of  the  element. 
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7.  Warping  of  the  beam  due  to  torsion  is  neglected. 

8.  The  constituitive  relations  for  the  plate  and  beam  material 
incorporate  the  von  Mises  yield  condition,  the  Prandtl-Reuss  flow  rule 
and  isotropic  work  hardening.   The  material  behavior  is  specified  by  a 
piecewise  linear  uniaxial  stress- strain  relationship.   The  plate  cross- 
section  is  divided  into  four  layers  for  the  stress  computations. 

9.  Offset  points  are  used  to  relate  beam  element  nodal  to  plate 
element  nodal  deflections. 

10.  An  initial  state  of  stress  feature  is  available  wherein  each 
element  will  be  in  a  uniform  state  of  stress  with  all  deflections  zero. 

11.  Permissible  loadings  are  distributed  (pressure)  loading  normal 
to  the  deformed  structure,  either  continuous  or  impulsive,  impulsive 
concentrated  forces  at  the  node  points,  and  an  automated  pressure  frag- 
mentation impulse  loading  from  the  NOL  BLAST  Code. 

12.  Account  is  taken  of  randomly  located  small  perforations  of  a 
plate  element  due  to  fragment  penetration  by  reducing  the  membrane  and 
bending  stiffnesses,  the  pressure  load,  and  the  mass. 

13.  The  timewise  numerical  integration  scheme  is  an  explicit  proce- 
dure . 

Ik.      The  lumped  mass  technique  is  used  to  develop  the  mass  matrix. 
Rotary  inertia  is  considered. 

15.  An  eigenvalue  routine  is  incorporated  to  automatically  deter- 
mine the  numerically  stable  time  increment. 

16.  The  program  has  a  restart  capability. 

17.  The  code  is  operational  on  the  CDC  6600  (FORTPAN  EXTENDED)  sys- 
tem.  A  structure  with  a  maximum  of  300  joints  and  400  members  can  be 
analyzed.   A  total  core  storage  of  225,000  words  OCTAL  (CDC  system)  or 
270,000  bytes  (IBM  320- 165  system)  is  required. 
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The  modifications  to  the  BR-1  Code  to  make  it  applicable  to  bullet 
and  fragment  penetration  into  fuel  tanks  can  be  divided  into  two  main 
categories:   (l)  incorporation  of  the  hydraulic  ram  pressure  loading, 
and  (2)  accounting  for  the  fluid- structure  interaction  phenomena.   Due 
to  the  difficulties  involved  in  treating  the  fully  coupled  problem,  the 
fluid- structure  interaction  phenomena  will  be  approximated  by  the  piston 
theory  described  in  II- A.   Thus,  since  the  BR-1  Code  does  not  account 
for  damping,  it  will  be  necessary  to  modify  the  code  by  adding  the 
damping  term  pew  to  the  equations  of  motion. 

In  the  piston  theory,  the  pressure  on  the  structure  is  p  +  pcv. 
where  p  is  the  pressure  in  the  incident  wave,  p  is  the  fluid  density, 
c  is  the  acoustic  velocity  in  the  fluid,  and  v.  is  the  fluid  velocity. 
Lundstrom  has  developed  a  computer  code  that  will  predict  the  pressure 
and  velocity  throughout  a  rectangular  volume  of  an  irrotational,  ideal 
compressible  fluid  due  to  a  penetrating  bullet.  The  boundaries  of  the 
fluid  may  be  either  free  or  rigid.  This  code  will  be  used  to  determine 
the  pressure  on  the  fuel  tank  as  input  data  to  the  BR-1  code. 
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D.   Analysis  of  Cavity  Motion 

In  order  to  gain  insight  to  forces  and  velocities  associated  with 
cavity  motion  two  problems  were  formulated  and  solved  by  Holm  (Ref. ,7). 
One  was  the  collapse  of  a  spherical  cavity,  an  event  which  occurs  when 
a  vacuum  exists  in  the  cavity.   The  other  was  the  oscillation  of  a 
spherical  cavity  with  entrapped  gases. 

Cavity  Collapse 

The  equations  of  motion  for  an  incompressible  fluid  or  density  p  can 

be  manipulated  into  the  form 

R 


fk  •  $ 

o 


b8-^ /{?-;=]  R2d 


where  R  is  the  cavity  radius.   The  pressure  within  the  cavity  is  p,  and 
at  infinity  p  .   At  time  t  =  0  the  radius  is  R  .   For  the  cavity  collapse 
problem  p  =  0,  and  equation  (l)  can  be  manipulated  into  the  following 
form 


t     -Ro  /  X^dX 


m^ 


X3 


(II-D-2) 


x 


where  X  =  R/R  .   Note  that  the  factor  multiplying  the  integral  has  the 
dimensions  of  time  and  defines  a  characteristic  time 

R 
o 

T  =    


'2  P«> 


3  p  (n-D-3) 


Two  examples  are  shown  in  Figures  1  and  2.   For  both  cases  the  collapse 
time  is  t  =  0.6  t.   In  terms  of  the  characteristic  time,  the  curve  in 
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either  figure  is  a  universal  curve  for  this  problem.  As  a  rule- of- thumb, 
cavity  collapse  is  about  60$  of  the  characteristic  time. 

Cavity  Oscillation 

When  p  /  0  in  Equation  (l),  there  occurs  oscillations  of  the  cavity. 
It  is  assumed  that  the  compression  and  expansion  of  the  trapped  gas  is 
adiabatic  and  lossless,  i.e.,  isentropic.   Consequently,  pressure  and 
cavity  radius  are  related  by 

P0  \R/  (II-D-10 

where  p  is  the  gas  pressure  when  R  =  R  .   The  ratio  of  heat  capacities 
C  /C  is  symbolized  by  y.      Introduction  of  Equation  (h)   into  Equation 
(l)  yields 


1       f  d2 


Z2  [z3  (a  +   1)  -  aZ3^  -2 


(H-D-5) 


where  a  new  parameter  a   has  been  introduced.   This  parameter  depends  on 
initial  ratio  of  pressures  within  the  cavity  and  at  infinity  and  is 
defined  by 

P^ 


a   = 


(Y"l)Poo  (II-D-6) 


In  Equation  (5)  Z  equals  R  /R  which  should  be  noted  is  the  reciprocal 
of  X  in  Equation  (2).   A  solution  of  Equation  (5)  is  shown  in  Fig.  3. 

The  complete  curve  cannot  be  obtained  from  a  single  integration.   When 

3  3v 

Z  (a+  l)  -  o;Z   -1  =  0,  the  integrand  approaches  infinity.   This 

occurs  when  a  maximum  or  minimum  in  radius  or  pressure  is  attained. 
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At  this  point  a  new  set  of  R  and  p  values  are  introduced  and  inte- 
^  o      o 

gration  proceeds. 

For  the  example  of  Fig.  3,  the  pressure  starts  at  p  /p  =6.67  and 
Z  =  1.   When  p/p  =  0.4  at  Z  =  .80  the  integrand  approaches  infinity. 
This  is  at  the  minimum  of  the  p/p  trace.   One  continues  at  this  point 
letting  p  /p  =0.4  and  a  new  Z  defined  as  1.000.   Integration  then 

O    00 

continues  to  get  the  portion  of  the  curve  from  minimum  p/p  to  maximum 

p/p  ,  i.e.,  the  part  of  the  curve  with  positive  slope.  For  the  case 

of  a   =  16.675?  one  cycle  of  oscillation  requires  a  time  of  1.35T. 

In  the  interpretation  of  Fig.  3  note  that  large  p/p  occurs  for 

small  R  or  large  Z.   Although  0.4  ^  p/p  ^  6.67  in  Fig.  3,  the  range  of 

radius  is  0.64  ^  R/R  <  1.24.   Large  variations  in  p/p  are  caused  by 
small  changes  in  R/R  due  to  3 7  power  dependence. 
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III.      HYDRAULIC   RAM  EXPERIMENTAL  FACILITIES 

A.   Introduction 

The  design,  fabrication,  and  calibration  of  a  ballistic  range  for 
aircraft  survivability  studies  was  accomplished  under  the  present  con- 
tract.  The  range  was  designed  specifically  with  the  hydraulic  ram 
problem  in  mind.   Of  paramount  importance  was  a  consistently  accurate 
method  for  determining  a  time  history  of  the  projectile  velocity  before 
projectile-tank  impact  and  during  hydraulic  ram  events  that  occur  after 
penetration.   Since  projectile  attitude  at  entry  is  important,  shadow- 
graph stations  were  required  along  the  flight  path  to  make  an  estimate 
of  projectile  attitude.   Two  shadowgraph  stations  were  also  required  to 
picture  the  hydraulic  ram  phenomenum  as  a  function  of  time  within  the 
tank.   The  shock  wave  system,  projectile  position,  and  cavity  shape  and 
size  could  be  determined  using  this  data.   The  electronics  associated 
with  the  ballistic  range  was  designed  to  give  stable,  repeatable  measure- 
ments at  very  small  time  intervals.   Calibration  of  the  present  range 
components  has  shown  that  consistently  accurate  data  is  obtained.   As 
a  result,  a  reliable  means  of  obtaining  projectile  energy  as  a  function 
of  time  during  hydraulic  ram  has  been  accomplished. 
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B.   Ballistic  Range  Components 

Basic  elements  of  the  ballistic  range  are  shown  in  Figure  III-B-1. 
The  rifle  mounting  system  is  composed  of  the  rifle  mount  and  stand.   The 
rifle  mount  was  adjustable  in  azimuth  (+  8  degrees)  and  elevation 
(+3  degrees  to  -  5  degrees)  for  ease  in  boresighting .   Rifles  of  22.2 
and  30  caliber  may  be  used.   Figure  III-B-2  shows  the  complete  rifle 
mounting  system  and  Figure  III-B-3  shows  the  rifle  mount  in  detail. 

The  shadowgraph  stations  are  composed  of  a  bullet  sensor,  time  delay 
unit,  spark  source,  collimating  lens,  and  a  shadow  box  with  reference 
grid  for  mounting  the  Polaroid  film  holder.   Figure  III-B-4  shows  a 
typical  shadowgraph  station.   The  bullet  sensor  is  a  chronograph  screen 
that  has  a  five  volt  D.C.  signal  shorted  to  ground  across  it.   The 
bullet  sensors  provide  start  and  stop  pulses  to  the  counters  and  delay 
circuits.   The  delay  units  generate  time  delayed  pulses  ranging  from 
190  lis  to  1700  u.s  that  trigger  spark  sources.   Figure  III-B-5  shows  the 
delay  unit  circuitry.   Figures  III-B-6  and  -7  shows  typical  shadowgraphs 
taken  along  the  bullets'  trajectory  before  impact.   The  complete  ballis- 
tic range  configuartion  is  shown  in  Figures  III-B-8  and  -9-   The  bullet 
catcher  is  shown  in  Figure  III-B-10. 
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Figure    IH-B-?        Ballistic   Ranffe   Rifle  -stern 
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Figure  lll-B-4    Ballistic  Range  Shadowgraph  Station 
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Figure  III-B-6 


Shadowgraph  of  a  63gr  Semi-Point  with  a 
Velocity  of  2941  f ps . 


Figure  III-B-7 


Shadowgraph  of  a  50gr  Spitzer  Used  in 
Range  Calibration  with  a  Velocity  of 
3080  fps. 
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Figure  III-B-9   Ballistic  Range  Configuration  (Up  Range  View) 
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C .   Experimental  Apparatus 

The  experimental  program  was  divided  into  three  phases.   The  first 
phase  consisted  of  measurements  of  the  projectile  energy  loss  during 
impact  with  the  wall  of  an  empty  tank.   Phase  II  measurements  were  taken 
with  a  water- filled  tank  but  with  the  projectile  entering  the  tank 
through  a  prepunched  hole.   Phase  III  measurements  were  taken  using 
solid  entry  walls  with  a  fluid- filled  tank.   A  comparison  of  these 
results  make  it  possible  to  partition  the  energy  losses  of  each  of  the 
various  hydraulic  ram  component  dissipative  processes.   Figure  III-C-1 
shows  the  test  apparatus  for  Phase  I.   A  test  plate  is  positioned  with 
two  bullet  sensors  mounted  directly  behind.   The  average  velocity  after 
plate  penetration  is  measured  and  the  interaction  phenomenum  may  be 
studied  using  a  shadowgraph. 

Phase  II  and  III  testing  uses  the  test  setup  shown  in  Figure  III-C-2, 
Two  shadowgraph  stations  provide  for  a  photographic  record  of  hydraulic 
ram  phenomenum  at  two  preset  times  after  projectile  impact. 
Figure  III-C-3  shows  the  plexiglass  side  wall  tank  with  its  removable 
entry  and  exit  walls.   The  entry  and  exit  walls  are  fastened  to  the 
1  inch  plexiglass  side  walls  by  studs  that  pass  through  steel  clamping 
frames  which  sandwich  the  test  plate.   The  tank  is  a  cube  with  18  inch 
inside  dimensions.   The  electronic  delay  and  spark  source  trigger  cir- 
cuits used  for  the  tank  shadowgraph  are  similar  to  those  shown  in 
Figure  III-B-5 •   Station  one  timing  starts  at  the  last  downrange  shadow- 
graph station  sensor  and  its  delay  circuit  range  is  variable  from  200  ^s 
to  165O  us.   Station  two  timing  starts  when  the  station  one  spark  source 
fires  and  another  delay  circuit  is  used  to  fire  the  station  two  spark 
source  from  3  M-s  to  20  [is   later. 
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Fig.    III-C-2     Fuel  Tank  Test  Setup 
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Fig.  III-C-3  Test  Tank  and  Stand 


IV.   TEST  RESULTS  TO  DATE 

A.   Penetration  Studies 

Phase  I  testing  was  conducted  to  gain  some  useful  knowledge  of  the 
fuel  tank  wall- projectile  penetration  characteristics.   Metal  plates  of 
7075- T6  aluminum  of  thicknesses  o.05,  0.09,  and  0.l6  inch  were  used. 
Projectiles  of  three  masses  and  three  nose  shapes  were  fired  into  these 
plates  at  varying  initial  velocities  using  a  .222  Remington  rifle.   The 
projectile  velocity  was  measured  before  and  after  the  plate  penetration. 
All  projectiles  studied  were  stable  before  impact  and  entered  the  test 
plate  perpendicular.   Several  shots  at  each  mass,  nose  shape,  and  velocity 
were  made  to  provide  statistical  data  in  order  to  determine  the  energy 
lost  by  the  projectile  during  penetration  of  the  entry  wall  without 
fluid  support.   The  energy  loss  thus  determined  will  be  the  same  as  the 
energy  lost  during  a  similar  event  in  hydraulic  ram  except  for  the 
effects  of  fluid  foundation  and  damping.   Phase  II  testing  will  yield 
the  amount  of  energy  initially  lost  to  the  fluid  with  entry  wall- fluid 
interactions  but  with  no  penetration  effects.   Phase  III  testing  will 
yield  the  total  initial  energy  loss  during  penetration.   Comparisons  of 
Phase  II  and  Phase  III  should  yield  the  energy  loss  during  entry  wall 
penetration  with  fluid  effects.   Comparison  of  these  results  with  Phase 
I  should  then  yield  the  difference  in  energy  loss  due  to  the  presence 
of  the  fluid  medium. 

To  date  the  Phase  I  testing  has  been  accomplished  (Ref.  8)  for  the 
0.09  inch  thickness  plate  only.  At  the  present  time  test  data  for  the 
other  two  wall  thicknesses  is  being  taken.  Plans  also  include  similar 
testing  using  30  caliber  projectiles  in  the  near  future. 

Figure  IV- A- 1  shows  a  typical  projectile  flowfield  after  plate 
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penetration.   Spallation  occurs  from  the  back  face  of  the  entry  wall 
producing  small  supersonic  fragments  that  travel  behind  the  projectile 
bow  shock.  A  nearly  spherical  strong  shock  wave  propagates  into  the 
tank  from  the  projectile  point  of  entry  during  the  initial  stage  of  plate 
penetration.   Shortly  after  bullet  passage  through  the  plate,  the  bullet 
bow  shock  wave  catches  up  to  the  original  spherical  wave.   These  waves 
coalesce  and  the  resulting  bullet  bow  shock  wave  has  a  reflex  shape. 

Previous  work  by  Forman  (Ref.  9)  derived  equations  for  the  minimum 
penetration  velocity  (ballistic  limit)  and  the  projectile  exit  velocity. 
The  ballistic  limit  is  given  by: 


o  -  y^i  d.6t  -  V3>  (iv-A-D 


V       15m 

where  G  is  the  material  shear  modulus,  r  is  the  radius  of  the  hole 
produced  in  the  plate  (assumed  to  be  the  same  as  the  radius  of  the  pro- 
jectile), t  the  wall  thickness,  and  m  is  the  projectile  mass.   For  impact 
velocities  greater  than  V  the  exit  velocity  V  approaches  the  impact 
velocity.   Assuming  a  constant  energy  loss  the  exit  velocity  may  be 
found  by: 


(IV- A- 2) 


where  V.  is  the  impact  velocity.   Equations  (l)  and  (2)  are  independent 
of  nose  shape  and  assumed  that  the  kinetic  energy  loss  is  independent  of 
impact  velocity  for  impact  velocities  greater  than  the  ballistic  limit. 
Figure  IV- A- 2  shows  the  results  of  firing  projectiles  of  constant  mass 
but  different  nose  shapes.   The  data  for  both  shapes  agree  closely  with 
the  theoretical  predictions.   This  indicates  that  nose  shape  has  a  minor 
role  in  determining  the  exit  velocity.   Figures  IV- A- 3  and  -k   show  the 
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result  of  tests  using  projectiles  of  different  mass  with  the  same  nose 
shape.   These  tests  again  showed  excellent  agreement  with  theory  over 
the  tested  range  of  impact  velocities.   The  greater  length  of  the  63  gr. 
projectiles  increase  the  effect  any  yaw  or  spin  axis  nutation  has  at 
plate  impact.   If  entry  is  not  perpendicular,  damage  is  more  severe  and 
energy  loss  due  to  penetration  is  greater.   A  plot  of  kinetic  energy  loss 
versus  initial  velocity  for  all  projectiles  is  shown  in  Figure  IV-A-5. 
The  theoretical  energy  loss  is  invariant  with  initial  velocity  and  is 
shown  as  a  reference.   The  energy  loss  associated  with  high  mass  project- 
iles is  typically  greater  than  that  of  the  lower  mass  projectile  for 
equal  initial  velocities  as  expected. 

Inspection  of  the  damage  to  the  entry  walls  shows  that  maximum 
damage  occurs  at  a  velocity  slightly  higher  than  the  ballistic  limit 
and  the  damage  size  approaches  the  diameter  of  the  bullet  as  the  impact 
velocity  increases  from  the  ballistic  limit. 
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Fig.  IV-A-2   r.xlt  Velocity  versus  Impact  Velocity  for  Two 
22.2  Caliber  Projectile  Shapes  Fired  Normal 
to  a  0.090  inch  Thick  7075-T6  Aluminium  plate 
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Fig.  IV -A -3 


Hxit  Velocity  versus  Impact  Velocity  for  a 
55gr  Semi-Pointed  22.2  Caliber  Projectile 
Fired  Normal  to  a  0.090  inch  Thick  7075-T6 
Aluminium  Plate. 
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Fig.  IV-A-U 


lixit  Velocity  versus  Impact  Velocity  for  a 
63gr  Semi-Pointed  22.2  Caliber  Projectile 
Fired  Normal  to  0.090  inch  Thick  7075-T6 
Aluminium  Plate. 
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B.   Shock  Phase  Studies 

Phase  II  and  III  studies  of  the  shock  phase  of  hydraulic  ram  to 
date  have  been  limited  to  low  energy  .22  caliber  shots.   An  attempt  to 
increase  kinetic  energy  level  by  using  higher  velocity  .222  ammunition 
resulted  in  catastropic  failure  of  the  test  tank  as  shown  in  Figure 
IV-B-1.   The  failure  of  this  tank  design  to  withstand  hydraulic  ram 
pressure  loadings  is  a  regrettably  graphic  example  of  the  magnitude  of 
the  survivable  fuel  tank  design  problem.   Inspection  of  the  wall  frac- 
ture shape  seems  to  indicate  that  hydraulic  ram  tensile  loads,  concen- 
trated locally  by  the  side  wall  stud  threads,  produced  failure.   It  can 
also  be  deduced  that  since  the  front  wall  sustained  the  major  damage, 
shock  phase  pressures  can  significantly  reduce  the  ability  of  fuel  tank 
front  walls  to  survive  the  later  phases  of  hydraulic  ram. 

In  spite  of  the  test  tank  destruction,  much  useful  information  has 
been  gathered  on  the  shock  phase  of  hydraulic  ram.   Dual  spark  shadow- 
graphs of  the  tank  were  taken  at  varying  times  after  projectile  impact. 
These  data  provided  information  on  the  characteristics  of  the  shock 
wave  pressure  pulse  and  the  rate  of  energy  loss  of  the  projectile. 

1.   Shadowgraphs  of  Fuel  Tank  Flow  Field 

Phase  II  testing  was  conducted  using  a  1  inch  diameter  pre- 
punched hole  in  a  .09  inch  thick  entry  wall.   Figure  IV-B-2  illustrates 
typical  results  of  the  impact  by  a  1-95  gram  lead  bullet  at  a  velocity 
of  0.37^-  kilometers  per  second  (1226  ft/ sec)  into  the  water- filled  test 
tank.   The  photographs  of  Figure  IV-B-2  show  that  a  hemispherical  wave 
front  is  produced  whose  center  is  the  point  of  impact.   It  is  implied 
from  this  fact  that  projectile  shape  and  forward  momentum  have  no 
significant  effect  on  the  wave  front  motion  during  the  early  stages  of 
hydraulic  ram.   Inspection  of  the  wave  thickness  also  shows  that  the 
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shock  front  strength  varies  along  the  wave.   The  greatest  strength 
occurs  in  the  direction  of  projectile  flight  and  the  least  strength 
occurs  where  the  shock  and  wall  meet. 

Phase  II  testing  involved  impacts  with  standard  factory  loaded 
.22LR  ammunition.   Average  velocities  of  the  resulting  shock  waves  were 
determined  from  shadowgraphs  similar  to  those  in  Figure  IV-B-2.   The 
shadowgraphs  were  used  to  determine  time  and  position  of  the  shock  wave 
as  it  propagated  through  the  fluid.   The  slopes  of  the  resulting  time- 
distance  plots  yielded  propagation  velocities  of  the  shock.   Figure  IV- 
B-3  in  such  a  plot  of  shock  front  velocity  for  a  projectile  with  an 
averaged  impact  kinetic  energy  of  135-2  Joules.   The  plot  shows  that  the 
shock  front  velocity  decays  to  the  acoustic  speed  in  water  (1.5  kilo- 
meters per  second  or  ^920  ft  per  second)  at  a  distance  of  about  11 
centimeters  (4.3  in). 

Pressures  generated  in  the  water  at  the  shock  front  corresponding 
to  the  measured  shock  front  velocities  were  determined  from  data  given 
in  Table  2.2  in  Ref.  10.   Figure  IV- B- h   shows  these  shock  front  pressures 
Extremely  high  pressures  are  generated  in  the  water  near  the  impact  point 
and  decay  rapidly  due  to  geometric  expansion  of  the  shock  front.   For 
example,  Figure  IV-B-4  shows  that  at  a  distance  of  2.2^  centimeters 
(0.88  in)  from  the  impact  point  the  resulting  shock  front  pressure  is 
3750  kilograms  per  square  centimeters  (53,335  lb/ ft  ).   As  with  any 
object  moving  through  a  fluid,  the  projectile  creates  a  dynamic  pressure 

in  the  water  due  to  the  stagnation  conditions  on  the  front  surface. 

h  2 

Pressures  of  710  kilograms  per  square  centimeter  (1.01  x  10  lb/ in  ) 

were  found  to  exist  in  the  vicinity  of  the  projectile  during  these 

studies.   Figure  IV-B-5  shows  the  evolution  of  the  pressure  field 

behind  the  shock  front  as  a  function  of  time.   When  the  projectile 


55 


first  penetrates  the  wall  a  strong  shock  is  created  and  decays  in 
strength  as  it  propagates  outward.   At  time  t  the  pressure  behind  the 
shock  greatly  exceeds  the  stagnation  pressure  on  the  projectile  nose. 
Comparing  pressures  at  t.  and  t  it  is  apparent  that  both  shock  front 
pressure  and  projectile  stagnation  pressure  decay;  however,  the  decay 
of  shock  front  pressure  is  much  more  rapid.   At  time  t  the  stagnation 
pressure  exceeds  the  shock  front  pressure.   One  could  define  the  time 
where  these  pressures  are  equal  as  a  characteristic  time  separating  the 
shock  and  cavity  phases  of  hydraulic  ram.   The  characteristic  time  for 
the  conditions  of  the  test  was  found  to  be  32  microseconds. 

Figure  IV-B-6  illustrates  the  later  stages  of  cavity  growth  and 
projectile  speed  decay.   The  air- filled  cavity  behind  the  projectile 
is  opaque  to  light  and  casts  a  sharp  shadow.   The  true  shape  of  the 
cavity  is  not  that  of  the  shadow  because  of  light  ray  refraction  through 
the  test  setup.   The  spark  shadow  of  the  cavity  is  wider  than  the  true 
cavity  size  which  could  be  found  using  X-ray  shadowgraph  techniques. 
The  projectile  shadow  is  not  visible  but  its  position  may  be  inferred 
by  the  position  of  the  bright  cusp  of  light  present  at  the  leading  edge 
of  the  cavity  shadow.   This  cusp  is  formed  by  light  bent  around  the 
projectile.   This  bending  is  caused  by  changes  in  fluid  refractive  index, 
It  is  expected  that  the  position  of  the  nose  of  the  projectile  corres- 
ponds to  the  tip  of  the  cusp.   Using  this  technique,  measurements  were 
made  of  the  progress  of  the  projectile  through  the  tank  after  impact. 

A  detailed  analysis  of  projectile  velocity  decay  is  extremely  diffi- 
cult because  of  the  many  variable  parameters.   Assumptions  must  be  made 
to  account  for  the  changes  in  projectile  shape,  drag  coefficient,  and 
attitude.   Because  of  these  difficulties  a  simplified  analysis  was  made. 
If  the  projectile  drag  is  given  by 
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the  assumption  of  constant  area,  drag  coefficient  equal  to  one  and 
constant  projectile  mass  results  in  the  following  projectile  velocity 
decay: 
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Figure  IV- B- 7  and  IV-B-8  show  a  comparison  of  experimental  data  with 
the  previous  analytical  results  for  velocity  and  energy  decay  as  a 
function  of  distance  traveled  by  the  projectile.   This  comparison  indi- 
cates that  in  spite  of  the  simplifying  assumptions,  the  agreement 
between  calculated  and  experimental  values  of  velocity  and  energy  after 
impact  was  good.   In  the  experimental  tests  the  projectile  underwent 
a  continual  deformation  during  deceleration  into  the  tank.   The  increase 
in  frontal  area  during  the  deceleration  was  approximated  by  using  an 
average  final  projectile  frontal  area  of  0.5  square  centimeters. 

Decay  of  the  projectile  kinetic  energy  can  be  directly  related  to 
the  amount  of  energy  being  transferred  to  the  fluid.   The  slope  of  the 
projectile  energy  curve  gives  the  rate  of  energy  transfer  to  the  fluid, 
assuming  no  projectile  deformation.   The  average  energy  transfer  rate 
for  the  projectile  tested  was  approximately  20.6  Joules  per  centimeter 
of  fluid  traversed. 

Phase  III  testing  was  accomplished  using  techniques  similar  to  those 
used  in  Phase  II.   The  test  entry  wall  was  however,  solid.   Figure  IV- 
B-9  illustrates  typical  results  of  the  impact  by  a  1.95  gram  lead 
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bullet  at  a  velocity  of  .37^  kilometers  per  second  (1226  ft/ sec).   The 
photographs  of  Figure  IV-B-9  show  the  obscuration  of  events  interior 
to  the  shock  wave.   The  shadowgraphs  show  a  strong  hemispherical  wave 
front  as  in  the  "pre-punched"  wall  case.   This  hemispherical  wave  is 
however,  interacting  with  another  wave  front  which  appears  as  a  straight 
line  oblique  to  the  wall.   This  wave  front  was  caused  by  the  propagation 
of  shear  and  dilatation  waves  moving  through  the  metal  wall  at  a  faster 
rate  than  the  waves  through  the  water.   This  wave  is  created  by  the 
penetration  stage  punching  action  of  the  projectile  and  its  rate  of 
propagation  was  found  to  be  proportional  to  the  bulk  sonic  velocity 
(dilatation  wave)  for  aluminum  (6.2  km/ sec.) 
2.   Deformation  of  Lead  Projectile 

Projectiles  were  recovered  after  each  experimental  test  and 
were  found  to  display  a  unique  pattern  of  deformation.   The  deformation 
phenomenon  was  very  consistent  for  projectiles  impacting  the  fluid 
through  pre-punched  plates.   Three  examples  of  projectiles  after  this 
type  of  impact  are  shown  in  Figure  IV-B-10.   The  extent  of  deformation 
was  directly  proportional  to  the  impact  velocity  and  all  projectiles 
exhibited  one  flat  side  and  an  undeformed  side.   This  type  of  deforma- 
tion is  the  result  of  the  dynamic  stability  of  the  projectile  as  it 
traverses  the  tank  fluid.   As  shown  in  Figure  IV-B-11  the  dynamic 
stability  of  the  projectile  is  maintained  through  the  equilibrium  of 
the  pressure  and  inertia.   The  figure  demonstrates  the  effect  of  an 
angular  displacement  of  the  projectile  from  the  equilibrium  position. 
The  center  of  pressure  moves  upward  on  a  slightly  curved  front  surface 
with  the  stagnation  point.   This  creates  a  restoring  moment  which  returns 
the  projectile  to  its  original  position.   Projectiles  which  penetrated 
a  solid  entry  wall  did  not  demonstrate  this  deformation  phenomenon. 
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Figure  IV-B-12  shows  these  projectiles  and  the  small  aluminum  cap  that 
was  punched  out  by  the  projectile  passage  through  the  entry  wall. 
Figure  IV-B-13  shows  typical  entry  wall  damage. 
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Figure   IV-B-1: 


Test  Tank  After   Impact  of 
High  Energy   Projectile 
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Figure  IV- B- 2: 


Shock  Waves  Resulting  from  a 
Prepunched  Hole  Impact 
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Figure   IV-B-3 


Shock  Front  Velocity  Generated  in  Water 
by   Projectile   Impact. 
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Figure   IV-B-A 


Shock  Front  Pressure  Generated  in 
Water  by  Projectile  Impact. 
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Figure  IV-B-5   Evolution  of  Pressure  in  the  Transition  from  Shock  to  Cavity  Phase 
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Figure  IV-B-6: 


Cavity  Growth  Resulting  from  a 
Prepunched  Hole  Impact 
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FIGURE  I2-B-7     PROJECTILE    VELOCITY    DECAY   RATIO 

AS    A    FUNCTION    OF   DISTANCE    INTO 
TANK 
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FIGURE  I2:-B-8     EXPERIMENTAL    AND    CALCULATED 

RATIOS    OF  PROJECTILE     ENERGY 
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Figure  IV-B-9 


Shock  Waves  Resulting  from  0.22  Caliber 
Camber  Projectile  Impact 
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Figure  IV-B-10 


Projectiles  After  Impact  Through 
Pre-punched  Entry  Wall 
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Figure   IV-B-ll 


Dynamic  Stability  of  a  Deformed 
Projectile  in  Water  After  Impact. 
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Figure  IV-B-12 


Projectiles  After  Impact  Through 
Solid  Entry  Wall. 
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Figure  IV-B-13   Solid  Entrv  Wall  After  Impact 
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V.   RECOMMENDATIONS  FOR  FURTHER  STUDY 

In  FY  7^,  Hydraulic  Ram  Studies  will  be  continued  in  the  Department 
of  Aeronautics.   The  following  items  are  recommended  as  an  outline  for 
this  future  work: 

1.  Shock  phase  shadowgraph  studies  should  be  extended  to  higher 
energy  and  caliber  projectiles.   Toward  this  end  an  experimental 
test  tank  capable  of  withstanding  the  damage  caused  by  these 
projectiles  has  been  designed.   Information  produced  by  the  above 
studies  should  be  correlated  with  analytical  predictions  of 
Yurkovitch. 

2.  The  test  tank  should  be  instrumented  to  measure  internal  pressures 
during  the  shock  and  cavity  phases.   This  data  will  allow  assess- 
ment of  the  Yurkovitch  predictions  during  shock  phase  and  - 
Lundstrom  predictions  during  the  cavity  phase . 

3.  The  test  tank  entry  wall  should  be  instrumented  with  strain 
gages  to  measure  wall  stresses  during  hydraulic  ram.   This  data 
will  allow  assessment  of  structural  interaction  prediction 
methods. 

k.      The  test  tank  entry  wall  should  be  instrumented  to  obtain  wall 
motion  during  hydraulic  ram.   A  fiber  optic  proximity  gage  is 
suggested  as  the  measuring  technique.   This  information  will 
allow  assessment  of  the  structural  interaction  prediction  method. 

5.   Continuing  studies  should  be  made  to  ascertain  the  validity  of  the 
combined  piston  theory  -  fluid  pressure  analytical- numerical 
approach  and  to  improve  the  methods  where  ever  possible. 
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